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Abstract 

Measurements of the kinematic distributions of J/xj) mesons produced in p— C, p— Ti and p— W collisions at 
— 41.6 GeV in the Feynman-a? region —0.34 < xf < 0.14 and for transverse momentum up to 
Pt = 5.4 GeV jc are presented. The xf and pr dependencies of the nuclear suppression parameter, a, are also 
given. The results arc based on 2.4 • 10 5 J/ip mesons in both the e + e and /i + fi decay channels. The data have 
been collected by the HERA-B experiment at the HERA proton ring of the DESY laboratory. The measurement 
explores the negative region of xf for the first time. The average value of a in the measured xf region is 
0.981 ± 0.015. The data suggest that the strong nuclear suppression of J/tp production previously observed at 
high xf turns into an enhancement at negative xf- 
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1 Introduction 

The DESY experiment HERA-B has measured in- 
clusive J/ip production in proton-carbon, proton- 
titanium and proton-tungsten collisions at a center- 
of-mass energy yfs = 41.6 GeV. The results are based 
on a sample of about 2.4 • 10 5 J/ip mesons recon- 
structed in both dilepton decay channels. A mea- 
surement of the atomic mass number dependence of 
J I ip production is derived from a comparison of the 
different samples. The atomic number dependence 
of inclusive particle production is often character- 
ized by a power law: u p a — & p n • A a where u p n 
is the proton-nucleon cross section and <t p a is cor- 
responding proton-nucleus cross section for a target 
of atomic mass number A. Previous measurements 
by E866 at Fcrmilab [1, 2] at y/s = 38.8 GeV and 
NA50 at CERN [3] at lower energy indicate that 
a <~ 0.94 — 0.95 at xf ~ and decreases to ~ 0.65 as 
xp approaches unity [2]. The results presented here 
provide a first measurement of nuclear effects in char- 
monium production extending into the negative part 
of the Feynman-x spectrum, —0.34 < xf < 0.14 1 . 

An understanding of the basic mechanisms respon- 
sible for the suppression of charmonium production 
in proton-nucleus collisions relative to proton-nucleon 
collisions is a prerequisite for the identification of pos- 
sible signals of new physics in high-energy heavy-ion 
data. Interpretations of the existing proton-nucleus 
data at positive xf rely on a delicate balance of 
several processes: nuclear absorption, shadowing of 
parton densities, energy loss, interactions with co- 
movers, hadronization of intrinsic cc components of 
the scattering nucleons, and so on. Ad hoc combi- 
nations of such elementary mechanisms, considered 
within various theoretical frameworks [4, 5, 6] are 
able to qualitatively reproduce the observed strong in- 
crease of J I ip suppression as xf approaches unity [2] . 

However the presently available data give little 
guidance in the largely unexplored negative-XF re- 
gion where other mechanisms such as formation-time 
effects can influence the effective nuclear path length 
of produced states [7, 8, 9, 10] and can potentially 
lead to a decidedly different behavior of a. In con- 
trast to J/ip production in the positive xp region, 
at negative xp the produced cc pair preferentially 
evolves into a charmonium state before leaving the 
nucleus and nuclear effects influencing the J/ip itself 
become important. Especially in this region, differ- 
ent models and approaches lead to contrasting pre- 
dictions, for example arising from differing assump- 



1 A slightly different range is quoted in previous HERA-B 
publications of Jjifi results due to minor differences in selection 
cuts. 



tions on energy loss of beam partons or produced cc 
pairs. By extending the current experimental knowl- 
edge of the nuclear dependence of J / ip production to- 
wards negative xf, the measurement described here 
provides new constraints for possible explanations of 
the observed nuclear modification pattern. The wide 
range of transverse momenta (up to 5.4 GeV/c) of the 
present measurement also permits a complementary 
measurement of the pr-broadening effect obtained at 
lower energies. 

The paper is divided into four main sections: an 
overview of the apparatus, trigger and data samples 
(Sect. 2), a description of the methods used for the 
selection and counting of J/ips (Sect. 3), the measure- 
ments of the kinematic distributions (Sect. 4) and the 
measurements of the nuclear dependence (Sect. 5), 
followed by concluding remarks (Sect. 6). 

2 Apparatus, data taking and 
Monte Carlo simulation 

HERA-B was a fixed-target experiment which studied 
particles produced in interactions of 920GcV/c pro- 
tons with wire targets positioned in the halo of the 
proton beam. The apparatus [11], shown in Fig. 1, 
was a forward spectrometer with acceptance ranging 
from 15 to 220mrad and from 15 to 160mrad in the 
bending (xz) and vertical (yz) planes, respectively. 
Because of this large acceptance and the fact that the 
J / ip decay tracks were measured before the muon de- 
tector (MUON) and the electromagnetic calorimeter 
(ECAL), HERA-B was the first fixed-target experi- 
ment with significant coverage in the region of nega- 
tive xf with an accessible range of xf € [—0.34, 0.14]. 

The target system [12] consisted of eight wires 
which were grouped into two stations separated by 
4 cm along the beam line. Each wire could be indi- 
vidually steered in the beam halo by a servo system in 
order to maintain a constant interaction rate. A total 
of five wires, differing in shape (round or rectangular), 
dimensions (between 50 jum and 500 /xm) and material 
(C, Ti and W) were used. Depending on running con- 
ditions and run-plan, cither a single wire or a pair of 
wires was active for any given data-taking run. The 
interaction rate was maintained in the range of 2 to 
6 MHz, depending on beam conditions and target. 

The vertex detector (VDS) [13] comprised eight 
planar stations of double-sided silicon micro-strip 
modules, seven of which were mounted in Roman pots 
built into the vacuum vessel and operated at a mini- 
mum distance of 10 mm from the beam. A track tra- 
versed typically three stations, yielding twelve mea- 
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Figure 1: Top view of the HERA-B detector. 



surement points in four stereo views. Vertex reso- 
lutions of 450 and 50 (im in the beam direction and 
in the transverse plane, respectively, were achieved. 
The eighth station was on a fixed mount immediately 
following the exit window of the main vertex system, 
2 m downstream from the target. 

The momenta of charged tracks were measured 
from their bending through a vertical magnetic field 
of integral 2.13 Tm. The main tracker was located 
between 2 to 13 m downstream of the target, with 
one station preceding the magnet, four stations im- 
mediately after the magnet (Pattern Chambers, PCs) 
and another two (Trigger Chambers, TCs) after the 
Ring Imaging Cherenkov detector (RICH). Each sta- 
tion contained an inner part [14] (made of micro-strip 
gas chambers and covering angles less than 20 mrad) 
which was not included in the trigger system and 
therefore does not play a role in the analysis pre- 
sented here. The outer part (OTR) [15] was com- 
posed of honeycomb drift chambers, with wire pitches 
of 5 mm closer to the beam pipe and 10 mm elsewhere. 
The final momentum resolution for muons was found 
to be a p /p[%] = (1.61 ± 0.02) + (0.0051 ± 0.0006) • 
p [GeV/c] [15]. 

The identification of the J/ij) in its dilepton de- 
cay modes as well as the first stage of the trigger 
system relied mainly on the signals provided by the 
ECAL [17] and MUON [18] systems. The ECAL was 
a sampling calorimeter using shashlik technology with 
Pb and W absorbers sandwiched between scintillator 
layers. It was divided into three sections (Inner, Mid- 
dle and Outer) with cell widths of 2.2, 5.5 and 11 cm, 
respectively, to roughly equalize occupancies. The in- 
ner section used W absorbers, while the middle and 
outer sections used Pb. The design was optimized 



for good electron/photon energy resolution and for 
electron-hadron discrimination. The final energy res- 
olution reached by the detector can be written in the 
form <t e /E = A/\/E®B (E measured in GeV), with 
A = 0.206, 0.118 and 0.108 and B = 0.012, 0.014 
and 0.014, for the Inner, Middle and Outer sections, 
respectively. The spatial resolution ranged from 1 to 
10 mm depending on the calorimeter section and on 
the energy of the particle [17]. 

The MUON system consisted of four tracking sta- 
tions interleaved with iron or concrete absorbers. As 
in the main tracker, two different technologies were 
used: gas pixel chambers in the innermost region and 
conventional tube chambers in the outer part. For 
the last two stations, not only the anode wires of the 
tubes but also the segmented cathodes were read out. 
Signals from the cathode pads were also given as in- 
puts to the trigger. 

The RICH [16] detector relied on a C4F10 radi- 
ator and was used extensively in other analyses for 
ir/K/p separation. It also played a small role in the 
present analysis as a means to reject backgrounds in 
the dilepton analysis, particularly from kaon decays. 

The trigger system selected both e + e~ and 
signatures and was organized into three levels: a pre- 
trigger [17, 19], a First Level Trigger (FLT [20]) and a 
software-based Second Level Trigger (SLT [21]). The 
pretrigger used signals from the ECAL and MUON 
detectors and required the presence of at least two 
reconstructed ECAL clusters with transverse energy 
above 1.1 GeV or the presence of two muon candi- 
dates, defined as coincidences of projective pads in 
the last two MUON layers. Starting from pretrigger 
seeds, the FLT attempted to find tracks in a subset of 
the OTR tracking layers and required that at least one 



of the seeds resulted in a reconstructed track. Start- 
ing again from the pretrigger seeds, the SLT searched 
for tracks inside regions-of-interest generated by the 
pretriggers using all OTR layers and continued the 
tracking through the VDS. Finally, at least two fully 
reconstructed tracks consistent with the hypothesis of 
a common vertex were required. Events passing the 
SLT were transferred to a computer farm which pro- 
vided full online reconstruction of a fraction of the 
events for data quality monitoring. The global trig- 
ger suppression factor, 5 x 10 4 , resulted in an event 
archival rate of about 100 Hz. 

A total of 160 million dilepton triggered events 
were recorded between October 2002 and February 
2003, together with an approximately fixed 10 Hz rate 
of minimum bias events which were used for monitor- 
ing and luminosity determination. The event samples 
were distributed between three target materials: car- 
bon (64%), tungsten (32%) and a small fraction with 
titanium (4%). 

A full Monte Carlo (MC) simulation is used to de- 
termine the triggering (except for FLT, see below), 
reconstruction and selection efficiencies. In view of 
the range of physics topics addressed by the exper- 
iment (pA inelastic interactions, meson decays and 
heavy flavor physics), the MC generator is built as a 
combination of two standard tools: Pythia 5.7 [22] 
for heavy flavor (b or c) quark production in pN in- 
teractions and subsequent hadronization and Fritiof 
7.02 [23] for light quark production, secondary inter- 
actions in detector materials and generic pA inelastic 
interactions. The production of the J/ip is simulated 
by first generating the basic hard process pN — > ccX 
and subsequent cc hadronization with Pythia and 
then giving the remaining energy and momentum (X) 
of the interaction to Fritiof for generation of fur- 
ther interactions inside the hit nucleus. The gener- 
ated particles are then given as input to the Geant 
3.21 based package [24] for full simulation of active 
(instrumented) and inactive (support structure) ele- 
ments of the detector and for the digitization of the 
electronic signals. 

To describe the kinematic characteristics of the 
produced J ftps as accurately as possible, an xf, Pt, 
and decay-angle dependent weight is assigned to each 
event and used in the subsequent analysis to force the 
simulated J/ip production and decay distributions to 
agree with the corresponding measured distributions 
for each target material. The weights are determined 
by an iterative procedure in which computed correc- 
tions are based on comparisons of MC event distribu- 
tions after reconstruction and selection cuts with the 
corresponding distributions from data. 



The FLT efficiency is derived from an efficiency 
map which is determined from the data itself. Since 
the SLT result is completely independent of the FLT, 
and since the FLT triggered on only one of the two 
lepton tracks, the efficiency map can be determined 
by a so-called tag-and-probe method. Using the ef- 
ficiency map, each event is assigned a weight which 
multiplies the kinematic weight discussed above. 

To accurately reproduce the actual working con- 
figuration at the time of data-taking and properly ac- 
count for time variations of working conditions, the 
full data taking period is divided into five calibra- 
tion periods of similar lengths, each matched by a 
corresponding simulation sample for which the effi- 
ciencies of the individual detector cells are evaluated 
and given as input to the MC. The MC samples are 
reconstructed and analyzed with the same methods 
and software packages used for the analysis of the 
real data. 



3 J/ip selection and counting 

The electron and muon candidates are selected with 
common track- and vertex-selection criteria while 
channel-specific methods are applied for lepton iden- 
tification and the treatment of the J/ip signal. 

All tracks passing the SLT are initially considered 
as lepton candidates. The track reconstruction proce- 
dure consists of finding straight segments in the VDS 
and PCs independently, matching them to each other 
and also to segments in the TCs. A full, iterative 
fit of found tracks is then performed. To reject in- 
correctly reconstructed or ghost tracks, loose cuts on 
the minimum number of hits in the VDS and OTR, 
as well as on the y 2 probability of the track fit are 
applied by the event reconstruction algorithms. The 
detector acceptance and trigger requirements effec- 
tively limit the momentum and transverse momen- 
tum ranges of lepton tracks to 5 < p < 200GeV/c 
and 0.7 < pr < 5.0GeV/c, respectively. For each 
event, among all possible pairs of oppositely charged 
lepton candidates consistent with a common vertex 
(x 2 probability greater than 1%), only the pair with 
the best particle identification (see below) for both 
leptons is accepted. 

3.1 Dimuon channel 

Since muons are the only particles having a significant 
probability of penetrating through the absorbers of 
the MUON detector, only minimal selection cuts are 
needed to obtain a clean sample (see also [25]). The 
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Figure 2: The (a) and e + e~ (b) invariant mass 

distributions for the full data sample in the J / ip mass 
region. The continuous lines represent the result of 
fits performed with the functions described in the 
text; the dashed lines are the fitted backgrounds. 



background of muons from pion and kaon decays is 
reduced by imposing tighter cuts on the quality of 
the track fit and on the matching of track segments 
in the VDS, OTR and MUON. Doing so rejects the 
typical "broken trajectories" produced by decays into 
leptons. Contamination from kaons is further reduced 
by discarding tracks with high values of the corre- 
sponding RICH likelihood. After all selection cuts, 
the background under the J/tp signal is reduced by a 
factor of 2.5 with respect to the triggered data with 
a loss of about 11% of the signal. 

Fig. 2(a) shows the resulting dimuon mass spec- 
trum together with the result of a fit to a sum of three 
functions [26], which model the J/ip and ip(2S) signals 
and the exponential background. The J/ip and ip(2S) 
signals are each modeled as a superposition of three 
Gaussians with a common mean which takes into ac- 
count track resolution and effects of Moliere scatter- 
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Table 1 : The numbers of reconstructed J/ips after all 
selection cuts in the dimuon and dielectron channels, 
and for different target materials. 



ing and a function representing the radiative tail due 
to the decay J/ip — > /i+^i~7 [26]. The background is 
described by an exponential of a second-order poly- 
nomial. The fitted position and width of the J/ip 
peak are 3.0930 ± 0.0002 GeV/c 2 and 40 ± 1 MeV/c 2 , 
respectively. 



3.2 Dielectron channel 

The J/ip selection in the dielectron channel is affected 
by major background contributions from charged 
hadrons which produce energetic ECAL clusters and 
by overlapping photon and charged-hadron energy de- 
posits in the ECAL. For this reason, the electron iden- 
tification requirements were the subject of careful op- 
timization studies which resulted in much more strin- 
gent selection cuts than for the muon sample. 

A cut on the transverse energy of the ECAL clus- 
ter (Et > 1.15 GcV) is applied in order to mask dif- 
ferent threshold cuts applied at the pretrigger level 
for the various acquisition periods. 

The reconstructed momentum vectors of electrons 
and positrons are corrected for energy loss from 
bremsstrahlung (BR) emission in the materials in 
front of the magnet. For each electron track, an at- 
tempt to identify an ECAL cluster due to a BR pho- 
ton is made by looking for an energy deposition in 
coincidence with the extrapolation of the associated 
VDS track segment to the ECAL. Any recovered en- 
ergy (about 18% of the initial electron energy on aver- 
age) is then added to the momentum measured by the 
tracking system. Since BR is a clear signature for an 
electron track, it is also exploited to obtain substan- 
tial background reduction which is essential for accu- 
rate J/ip counting. The baseline results of the anal- 
ysis are obtained by requiring that at least one lep- 
ton of a decaying J/ip has an associated BR cluster. 
This requirement reduces the signal by about 30% and 
suppresses the background by more than a factor of 
two. Alternative requirements (no BR requirement, 
only one, or both electrons emitting BR) lead to very 
different background shapes and amounts. The dif- 
ferences are exploited for systematic studies on the 
stability of signal counting and on the correctness of 




Figure 3: Distribution of the number of J/ips reconstructed in the carbon target sample (points) and in the 
corresponding MC data (histograms, arbitrarily re- normalized) as a function of pr and xp. (a) and (b): muon 
channel, (c) and (d): electron channel. 



the MC simulation. 

Adjustments to the measured momenta of electron 
tracks were applied to compensate for differences in 
multiple scattering between electrons and muons since 
the track-fitter had been calibrated for muons. For 
this purpose, a correction map, determined from the 
shift of the J/ip peak position in different kinematic 
regions was used. 

Additional selection cuts are applied to further 
improve the significance of the dielectron signal. A 
particularly discriminating variable is the ratio E/p, 
where E is the energy of an ECAL cluster and p is 
the momentum of the associated track. The E/p dis- 
tribution for electrons has a Gaussian shape with a 
mean value close to 1 and width varying between 6.4% 
and 7.4% depending on calorimeter sector. Values of 
E/p much lower than 1 are mainly due to particles, 
mostly hadrons, which release only part of their en- 
ergy in the calorimeter. Further selection variables 



used in the analysis are the distances Ax and Ay - 
along the x and y directions - between the recon- 
structed cluster and the track position extrapolated 
to the ECAL. The Ax and Ay distributions for elec- 
trons are, apart from a small tail, well described by 
Gaussians centered at zero with widths between 0.2 
and 1.0 cm depending on calorimeter sector. Cuts 
on these quantities lead to a significant reduction of 
the contamination from hadrons and random cluster- 
track matches which are characterized by significantly 
wider distributions. The selection of the candidate 
electron-positron pairs is further refined by putting 
an upper bound on the distance of closest approach 
(Ab) between the two accepted tracks near the vertex. 

All the requirements described above have been 
simultaneously optimized by maximizing the signifi- 
cance S/^JS + B of the J/ip signal (S) - taken from 
the MC (scaled to the number of J/ip in the data) 
- with respect to the background (B) - evaluated 



from the data. The optimal ranges for the different 
cut variables depend on the number (one or two) of 
BR clusters associated to the electron-positron pair. 
When both electron and positron have a BR clus- 
ter correlated to the track, where the cluster posi- 
tion was determined by hierarchical clustering [17], 
the event is already rather cleanly reconstructed and 
only one additional request (for each lepton candi- 
date), (E/p — l)/a E / p > —3.6, is applied. When 
only one of the two possible BR clusters is found, 
the accepted ranges are determined for each lepton 
as -3.6 < (E/p - l)/a E/p < 5.4, \Ax\/a Ax < 7.0, 
|Ay|/fAj/ < 3.3 and Ab < 370 yum, respectively. 

The combined selection cuts increase the S/B ra- 
tio of the J/tp by about a factor of 10 with respect 
to triggered events, and have an overall efficiency of 
(45 ±4)% as evaluated using the data and verified 
with the simulation the rather large uncertainty is 
due to the difficulty of counting the signal when no 
cuts are applied. As can be seen in Fig. 2(b), the sig- 
nificance of the optimized J/tp signal, although less 
than that of the muon channel shown in Fig. 2(a), is 
nonetheless such that the electron sample significantly 
enhances the statistical precision of the final results. 
The method adopted for counting the signal uses a 
Gaussian shape for the right part of the peak and 
a Breit-Wigner shape for the left part to account for 
the sizable asymmetry of the signal caused by missing 
BR energy and the contribution of the radiative decay 
J/tp — > e + e~7. The background is parametrized with 
a Gaussian at lower mass values and an exponential 
at higher mass, joined together such that the result- 
ing curve is smooth. The position and width of the 
J/tp peak as determined from the fit are 3.110±0.001 
GeV/c 2 and 72 ± lMeV/c 2 , respectively. 




p T (GeV/c) 

Figure 4: Inclusive pr distributions of J/tp mesons for 
three target materials with arbitrary normalizations. 
The error bars represent the combination of statisti- 
cal and systematic uncertainties. The interpolating 
lines are the results of a simultaneous fit of the three 
Pt distributions to Eq. 2 performed with the method 
described in Sec. 4.2. 



The yields of selected J/tp candidates for the two 
decay channels and for each target material are listed 
in Table 1. Fig. 3 shows (for the carbon data) a com- 
parison between data and MC of the distributions 
of reconstructed J ftps as a function of the kinematic 
variables pt and x p . 



4 Kinematic distributions 
4.1 Results 

The present analysis adopts the degrees of freedom 
Pt, xf and $ (azimuthal production angle) for the 
description of the J/tp production kinematics. Single- 
variable distributions are obtained according to the 



formula (here written e.g. for xf) 

dNj„ AN^xf)_ 

— j \ X F) - f \ T , (1) 

dx F £j/4,(xf) Ax F 

where AN r f^(xp) is the fraction of J ftps recon- 
structed in a given xp interval of width Axf and 
ej/^(xp) is the global (trigger, reconstruction and se- 
lection) efficiency for that interval integrated over all 
other kinematic variables (including the J/tp decay 
degrees of freedom) using the tuned MC. In each case, 
the signal is also integrated over all other kinematic 
variables. All distributions are normalized to unit 
area 2 . 

The final pt and xp distributions for the three 

2 The absolute J/ip yield in proton- nucleus collisions at 
920GeV/c was the subject of a measurement performed using 
minimum bias data [27] . 
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Table 2: J/ifj pt distributions (dN/dpr, normalized to their integrals over the measured range) for three 
target materials with statistical and systematic uncertainties. 



materials are shown in Figs. 4 and 5. The error bars 
include statistical and systematic uncertainties added 
in quadrature. The corresponding numbers can be 
found in Tables 2 and 3. The final distributions and 
systematic uncertainties are evaluated by fixing the 
input parameters and assumptions of the analysis to 
a variety of values within their range of uncertainties 
(see list below) and carrying through the full analysis. 
The kinematic distributions are finally obtained by 
averaging over decay channels for each set of input 
parameters and assumptions. The central value for 
each bin is the mid-point of the distribution of values 
thus obtained and the systematic uncertainty is the 
maximum spread of the obtained values divided by 
\/l2. The stability tests arc described in the following 
list. 

• The impact of selection and optimization re- 
quirements is evaluated by changing the cuts on 
momentum and transverse momentum of muon 
and electron candidates with respect to the in- 
trinsic thresholds of the trigger selection, and 



by scanning systematically the values of all cut 
variables used for the optimization of the dimuon 
and dielectron signals (including, for the latter, 
different BR requirements). 

• The uncertainty associated to the signal count- 
ing method has been estimated from the varia- 
tion of the results obtained with the adoption of 
modified background and signal functions. Spe- 
cial attention is given to the background eval- 
uation of the dielectron channel: as an alterna- 
tive to the fit with an assumed background func- 
tion, a background shape constructed by mixing 
real events (combining each track with one of 
opposite-charge from a different event) has been 
used in an unbinned maximum-likelihood fit of 
the invariant mass spectrum. A further cross- 
check is represented by the comparison between 
the efficiency-corrected J/ip yield obtained with 
different BR requirements (and therefore differ- 
ent background shapes). When the BR require- 
ment is removed, the evaluation of the number 
of J/ips becomes less stable due to increased 
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Table 3: J/ip xp distributions (dN/dxF, normalized to their integrals over the measured range) for the three 
target materials with statistical and systematic uncertainties. 



background, but the variation of the efficiency- 
corrected yield with respect to the standard se- 
lection is estimated to be lower than 5%. 

• The systematic uncertainties also account for the 
stability of the results when specific acquisition 
periods and conditions are selected. A large sub- 
sample of the collected events was produced on 
two target wires of different materials operated 
simultaneously. The comparison of these data 
with those acquired with a single target pro- 
vides an indication of the extent to which the 
measurements are affected by variations of the 
experimental conditions. 

• The results are sensitive to the shape of the J/ip 
decay angular distribution assumed in the MC 
generator. The hypothesis - made by previous 
experiments - that the J / ip is produced in an un- 
polarized state is not supported by the HERA-B 
data [28] . There is, moreover, an indication that 
the polarization increases in magnitude with de- 
creasing pt , while no significant x p dependence 
is found. Since a longitudinally polarized J/ip is 
detected more efficiently due to the lower prob- 
ability that its decay leptons escape detection 
by passing through the uninstrumented region 



near the beam, the kinematic dependence of the 
polarization assumed in the MC influences the 
shape of the efficiency-corrected pp and/or xf 
spectra. The systematic stability tests therefore 
include a variety of different assumptions for po- 
larization (including longitudinal, py-dependent 
polarization also considering the possibility of 
an yl-dcpendent polarization - and the absence 
of polarization). 

The distribution of the azimuthal production an- 
gle $ has been evaluated as a systematic check of the 
uniformity of the MC description of the geometrical 
acceptance. Fig. 6 shows the result obtained when 
combining the full data reconstructed in both decay 
channels: the points are consistent, within the statis- 
tical uncertainties, with the expected flat distribution. 

The dimuon and dielectron results (which, as dis- 
cussed below, are found to be compatible) are aver- 
aged such that correlations in their systematic un- 
certainties are taken into account. Such correlations 
have been estimated by maintaining (when possible) a 
parallelism among the two channels when evaluating 
the effect of each single systematic test. The assumed 
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Table 4: Parameter values obtained from the fit of the kinematic distributions of each of the three target 
samples to the functions described in the text (Eq.s 2 and 3). The first of the given uncertainty ranges is 
statistical and the second is systematic. 
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Figure 5: Inclusive xf distributions of J/ip mesons 
for the three target materials with arbitrary normal- 
izations. The error bars represent the combination of 
statistical and systematic uncertainties. The interpo- 
lating lines are the results of a simultaneous fit of the 
three xf distributions to Eq. 3 performed with the 
method described in Sec. 4.2. 
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Figure 6: Distribution of the azimuthal production 
angle of the J ftps fitted to a constant. The error bars 
include only statistical contributions. 



polarization hypothesis is found to be the dominat- 
ing source of uncertainty in the final results - espe- 
cially for the lower part of the pr distribution - as 
well as the most important cause of correlation be- 
tween the two analyses. Relative to this uncertainty, 
the signal selection cuts are in general responsible for 
negligible systematic variations, except for the most 
positive part of the x F spectrum, where acceptance 
corrections increase dramatically due to the low-angle 
detector acceptance cut-off near the beam. 

The results are well represented - and can there- 
fore be described - by the following interpolating 
functions which are further motivated in the following 
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Figure 7: The (pj) of the produced J/tps func- 
tion of A 1//3 — 1 (sec Eq. 4). The results of the present 
analysis (black filled circles with total and statisti- 
cal uncertainties) are compared to previous measure- 
ments performed with different beam energies at Fer- 
milab [1] and at the SPS [3]. The data are fitted with 
linear functions. 
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The parameters (pf.), /?, w XF (width at half maxi- 
mum), Axf (shift of the center of the distribution 
with respect to = 0) and 7 are left free in the fit 
of the distributions. The resulting values are listed 
in Table 4 for each target material. The systematic 
uncertainties of the parameters have been determined 
from the maximum variation (divided by \/l2) of the 
results obtained by re-fitting the distributions after 
each of the stability tests described above. The nor- 
malized x 2 obtained from the fits of the combined 
dimuon-diclcctron results are, respectively, 1.5, 0.8 
and 2.1 for the pt distributions of carbon, titanium 
and tungsten, and 1.1, 1.5 and 1.2 for the three Xp dis- 
tributions (only the statistical uncertainties are taken 
into account). 

Also included in Table 4 are the results obtained 
separately in the dimuon and dielectron channels with 
the respective systematic uncertainties. The good 
agreement between the results of the two analyses 
within the statistical uncertainties confirms that the 
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Figure 8: Energy dependence of the parameters 
(p^) pp (a) and p (b) as defined in Eq. 4 describing the 
pr-broadening of the J ftp. The points represent the 
results obtained by HERA-B (black filled circles, with 
total and statistical uncertainties), at Fermilab [1] 
and at the SPS [3]. The data are fitted with linear 
functions. 



channel-specific issues of particle identification and 
counting are not responsible for large systematic vari- 
ations in the shapes of the distributions. 



4.2 Parameterization and interpreta- 
tion of kinematic distributions 

Among the parameters adopted for the description 
of the data, the width of the pr distribution ((p^)), 
the position of the maximum of the xf distribution 
(Axf) and, possibly but less significantly, its width 
( w x F ) show a trend with the mass number A. 

It is well known [29] that the average p\ of parti- 
cles produced in nuclear collisions increases with the 
mass of the target nucleus. This trend is confirmed by 
the HERA-B data with high significance. The u Pt~ 
broadening" effect is commonly explained as a conse- 
quence of multiple elastic scattering of the incoming 
beam parton in the surrounding nucleus before the 
hard scattering process takes place. The measured 
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72/60 



0.00044 ± 
0.0017 ± 
0.00066 ± 
0.029 ± 



0.00022 
0.0046 
0.00092 
0.014 



Xf distribution, r = 



w>. 



X 2 /NDoF 



-0.00295 ± 0.00035 ± 0.00018 
0.1478 ± 0.0014 ± 0.0028 

(fixed) 
1.725 ± 0.0028 ±0.0015 
91.5/61 



Table 5: Results of the global fits of px and xf distri- 
butions described in the text. The first of the given 
uncertainty ranges is statistical and the second is sys- 
tematic. 



Figure 9: The width (a) and shift (b) of the J/ip xp 
distributions as a function of A 1 / 3 — 1 (see Eqs. 5 and 
6). The double error bars represent total and statis- 
tical uncertainties. The data are fitted with linear 
functions. 



increase of (p T ) with A is shown in Fig. 7 together 
with the results of experiments at lower energies. The 
variable of the abscissa, A 1 ! 3 — 1, is approximately 
proportional to the radius of the target nucleus i.e. 
to the average path length of the parton inside the 
nucleus with the shift of — 1 such that the magnitude 
of the effect is measured with respect to A = 1. All 
measurements are actually consistent with the param- 
eterization 

{p T )^{p T Y P + p{A 1 l 3 -l), (4) 

to which they are fitted in the plot. As summarized in 
Fig. 8(a), the results for the average p T extrapolated 
to proton-nucleon interactions ((p T ) pp ) are compati- 
ble with a linear growth with the square of the center- 
of-mass production energy s. On the other hand, p is 
approximately independent of center-of-mass energy 
as can be seen in Fig. 8(b). 

Furthermore, HERA-B observes a difference in 
shape between the xf distributions of the J/tp for 
different target nuclei consisting of an increasing dis- 
placement of the center of the distribution towards 
negative values. As shown in Fig. 9, J/ips are pro- 
duced in tungsten with an xp distribution which 



has equal or slightly greater width with respect to 
those produced in carbon and tends to be asymmet- 
rically centered at a lower value. This behavior is 
also supported by a fit of the E789 gold data [1] 
(-0.035 < x F < 0.135, E b = 800 GeV) with Eq. 3. 
The fitted width of 0.11 ±0.01 is lower than our value 
for tungsten suggesting not only that the maximum is 
shifted but that the shape becomes asymmetric. As a 
possible interpretation, the effect may be attributed 
to the energy loss undergone by the incident parton 
and/or the produced state in their path through the 
nucleus, causing a reduction of the average xf of the 
J/ip and a possible additional smearing of the mo- 
mentum distribution. This hypothesis motivates the 
choice of representing the data also in this case as a 
function of A 1 / 3 — 1. The points in Fig. 9 (a) and b) 
are fitted respectively with: 

w XF = wZ+t(A^ 3 -1), (5) 
Ax F = k {A 1 ' 3 - 1), (6) 

where k, vjPP, and r are free parameters. 

To obtain the best description of the dependence 
of the pp and xp spectra on the target nucleus, a 
simultaneous fit of the three (C, Ti and W) distribu- 
tions according to the functions given in Eq. 2 and 
3 has been done. According to the hypothesis that 
energy loss is responsible for the observed nuclear de- 
pendence of the shape of the kinematic distributions, 
the fit has been constrained by imposing the relations 



from Eqs. 4, 5 and 6, and, moreover, 

(3 = f3 pp + p'iA 1 / 3 - 1), (7) 
7 = 7 PP (independent of A), (8) 

where [3 PP , 0' and j pp are additional parameters of 
the fit. In Table 5 the results of this procedure are 
summarized. The fit of the xf distributions has been 
performed in two variants, with the parameter t left 
free or fixed to zero - therefore assuming in the lat- 
ter case that w XF is independent of A. The resulting 
best-fit curves (with r left free) are the interpolating 
lines plotted in Figs. 4 and 5. The fit results indicate 
a significant nuclear dependence not only of the pr 
distribution (parameter p), but also of the Xf distri- 
bution: there is a significance of la for k =/= when 
r is fixed to zero, which changes to 4 and 3er, re- 
spectively, for k and r (with a strong anti-correlation 
between the two) when both are left free. 

5 Nuclear dependence of J/ip 
production 

The Glauber Model [30] suggests that the dependence 
of the J/tp production cross section on atomic mass 
number (A) can be approximated by a power law: 

&pA — &pN ■ A a , (9) 

where <j p n is the proton-nucleon cross section and 
a, the "suppression" parameter, characterizes the nu- 
clear dependence. Pure hard scattering in the absence 
of any nuclear effects would correspond to a equal to 
unity. A suppression of J/tp production would lead 
to a < 1 while an enhancement (anti-screening effect) 
would be signaled by a > 1. Usually, a is described 
and measured as a function of xp and px (see for 
example [2, 3, 31]). 

Eq. 9 is generally used to describe data and pre- 
dictions independently of particular mechanisms of 
nuclear modification. In general, however, a may de- 
pend on A and thus depend on the targets used to 
make the measurement. For the measurement pre- 
sented here, a is evaluated by comparing the J/tp 
yields from two different targets: carbon and tung- 
sten 3 as a function of xf and pr- 



3 The titanium sample is not used for this analysis since 
it is too small to have a significant impact on the statistical 
precision of the result and would have required a significantly 
more complex analysis procedure. 



5.1 The a measurement 

Using Eq. 9, the nuclear suppression parameter a can 
be extracted from a measurement of the ratio of cross- 
sections for carbon (C) and tungsten (W) targets: 



GpC 




The measurement of the cross section ratio re- 
quires a measurement of the ratio of the integrated 
luminosities of the carbon and tungsten target sam- 
ples. For the HERA-B setup, this can be done using 
data samples where two different targets are oper- 
ated simultaneously (double-target runs) since most 
of the systematic uncertainties cancel and an abso- 
lute luminosity measurement can be avoided. On the 
other hand, for studies of the dependence of a on the 
kinematic variables, greater statistical precision can 
be obtained by also using the single-target runs. The 
HERA-B measurement of a thus consists of two sub- 
measurements: a measurement of the average value 
of a, (a), over the full visible kinematic range and 
a measurement of a — {a) as a function of xf and 
Pt- The shape distributions are then corrected us- 
ing (a) to produce an absolute measurement of the 
distribution of a over the measured range. 

More specifically, (a) is evaluated using double 
target runs based on the formula: 

_ M^) 1 { (Nw £c ec\ 

ln(^) ln(^)' U {n c 'C w 'e w J ' 

(11) 

where Nx (X=C,W) denotes the total number of re- 
constructed J/tp mesons originating from the corre- 
sponding target wire, Cx is the luminosity, ex is the 
overall detection efficiency (see Sec. 2) and the event 
yields are derived as discussed in Sec. 3. The measure- 
ment of the luminosity ratios is described in Sec. 5.2. 

The dependence of a on and pt is obtained 
from the full carbon and tungsten target data samples 
(double- and single-target runs). The full carbon sam- 
ple is roughly twice the size of the double-target sub- 
sample while the full tungsten sample is 10% larger 
than the double-target subsample. The shape of the 
differential distributions are given by (here written 
e.g. for xf)'- 

J_.^ (12) 
<jj/ip axF 

where (Jj/^p = <r(pA — > J/tp + X) is the total visible 
J/tp cross section. 

The measurement of nuclear effects can then be 
derived from the distributions of Eq. 12 using (here 



written e. g. for xj?): 



oi(xf) = — 7 r- • In 



ln (4t) 



1 _ daw 

<7\Y dxp 

1 da~ 

ac dx F 



+ (a) . (13) 



5.2 Luminosity ratios 

The luminosity acquired on target X, where X is ei- 
ther C (carbon) or W (tungsten), can be expressed 
as: 

N x A x • N BX 



T inel 
7 X 



where iVx is the total number of inelastic interactions 
occurring on target X during the measurement, a^ cl 
is the total inelastic cross section, N BX is the corre- 
sponding total number of filled bunch crossings (BX) 
and Ax is the average number of interactions per filled 
BX. The total cross sections a^ el for each target ma- 
terial together with some details on this topic can 
be found in [32]. The luminosity ratio Rc needed in 
Eq. 11 is then given by: 



T inel 
W 



Ac 



R £ - — - ■ — ■ 



-w 



(15) 



Assuming the interaction probability on target X 
follows a Poisson distribution, can be calculated 
from the observed number of events with at least one 
interaction (N(> l) obs ) using: 



1 



Ariel 



In 1 



obs 



(16) 



where e™ el is the probability to observe a single in- 
teraction. 

The determination of Ax relies on random-trigger 
events which were accumulated together with the 
dilepton-trigger events used for J/tjj counting. Five 
methods which differ by the event characteristics used 
to define the presence of an interaction are used to 
count events. All methods rely on tracks found in 
the vertex detector. To maintain high efficiency, the 
requirements imposed are minimal but sufficient to 
also keep the probability of incorrect target wire as- 
signment at a low level. The methods are based on 
the following five criteria: 

For all events, 

1 . > 1 primary vertex on wire X where the pri- 
mary vertex is formed from tracks measured 
both in the VDS and OTR ("long tracks"), 



2. > 2 tracks (including long tracks and tracks 
seen only in the VDS) with impact parame- 
ter < 3ct of wire X and > 5<r from the other 
wire, where a is the impact parameter mea- 
surement uncertainty, 

and, using the subset of events with no vertex 
found on the other wire, 

3. > 1 primary vertex on wire X using all 
tracks, 

4. > 1 primary vertex on wire X using long 
tracks only, 

5. > 2 long tracks within < 3cx from wire X. 

Both as a cross check and as an estimate of the 
systematic uncertainty on the efficiency, all counting 
methods are checked in parallel. For the final lumi- 
nosity ratio determination, the average of the five de- 
terminations is used and the rms spread of the five 
is factored into the systematic uncertainty. Further- 
more, Eq. 16 assumes that the interaction probabili- 
ties for each wire follow a Poisson distribution. How- 
ever, the individual bunch fillings are often uneven 
and the interaction rate varies in time by typically 
20%. To quantify this influence, an alternative lumi- 
nosity calculation is performed in which the detailed 
bunch filling structure and the interaction rate distri- 
bution on each wire are taken into account. The dif- 
ferences between the resulting ratios and those com- 
puted directly from Eq. 15 are negligible compared to 
other systematic uncertainties. 

To minimize the dependence of the efficiency es- 
timate on MC, the efficiency of each of the above 
methods is calibrated by comparing the luminosity 
estimate found using it with that found by the meth- 
ods described in [32]. These latter methods rely on 
very simple criteria to identify events with interac- 
tions, such as a minimal number of hits in the RICH 
detector (typically twenty, compared to thirty hits ex- 
pected for a fully accepted fast charged particle) or 
a small (1 GeV) energy deposit in the ECAL and are 
estimated to be sensitive to roughly 95% of the total 
non-diffractive cross section. 

Each target of a two-wire configuration is cali- 
brated separately using single-wire data runs taken 
nearby in time to the run being calibrated. A "ghost" 
wire is introduced at the location of the other wire 
of the configuration. Thus the MC is not relied on 
to model tracking in the VDS or vertex finding, but 
only to estimate the efficiency of simpler and more 
robust event counting techniques described in [32]. 
Overall efficiencies in the range of 60 - 80 % are 
found, depending on method and wire. The efficiency- 
calibration method based on real data is also used 
to evaluate the probability that interactions are as- 



signed to the wrong wire. This probability is method 
and configuration dependent and is never more than 
0.4%. 

The average relative systematic uncertainties on 
the luminosity ratios due to interaction counting and 
MC calibration [32] are 1.3% and 3.2%, respectively, 
giving an overall scale uncertainty of 3.4% on Re- 
Depending on wire configuration between 0.6 and 1.2 
million events were used for the determination of Ac 
and Aw, thus the statistical uncertainty on the lumi- 
nosity ratio is negligible. 



a 1.3: 
1.25 = 

1.2: 

1.15 = 
1.1 : 
1.05! 

1 : 
0.95 = 



HERA-B 2008 e + ]i 
E866 2000 



JO 



0.9 Mr^ 



0.85^ 



5 6 
Pt [GeV/c] 



5.3 Results 

Based on Eq. 1 1 , an average suppression value of 



Figure 10: The nuclear suppression parameter a as a 
function of pp measured by HERA-B (filled triangles, 
plotted with total and statistical uncertainties) and 
by E866 [2] (empty squares). 



(a) = 0.981 ± 0.004 stat . ± 0.016 sys . (17) 

in the visible range of xp is obtained. As explained 
in Sect. 2, the target system of HERA-B consisted 
of eight different wires grouped in two stations. The 
data were taken with four different two-wire configu- 
rations which were analyzed separately and averaged 
to obtain the final value. Also, the electron and the 
muon decay channels represent two statistically in- 
dependent measurements. The average value is cal- 
culated as a weighted mean with weights being the 
squared quadratic sum of statistical and luminosity 
ratio uncertainties. The luminosity ratios are deter- 
mined for each wire configuration separately as dis- 
cussed in Sec. 5.2 with a contribution to the sys- 
tematic uncertainty on (a) of 3.4 %/ h\(Aw/Ac) = 
1.24 %. A systematic effect of 1.1% due to time vari- 
atons of detector performance and imprecisions in de- 
tector or trigger simulations was estimated from the 
variations of (a) among the four samples and two de- 
cay channels. The total systematic uncertainty on 
(a) is thus 1.66 %. The statistical precision of the 
(a) measurement contributes an uncorrelated uncer- 
tainty 0.4 %. 

The values of a for individual x p and pp bins are 
given in Table 6 and shown in Fig. 10 for pp and 
in Fig. 11 for xp. They will be further discussed in 
Sec. 5.4. The values and systematic uncertainty es- 
timates were derived according to the procedure de- 
scribed in Sec. 4.1. The error bars on the figures show 
both statistical and total contributions. The system- 
atic uncertainties in the estimate of (a) are largely 
uncorrelated with those from the a — (a) measure- 
ment. The final systematic uncertainty estimate is 
the quadratic sum of the two. The systematic uncer- 
tainty is substantially correlated from bin to bin. 



5.4 Discussion 

The results of the a measurement as functions of pp 
and xp are given in Table 6. The distributions of 
the a(pp) an d a(xp) are presented in Figs. 10 and 
1 1 where they are compared with measurements per- 
formed by other fixed target experiments: E866 [2] 
{E p = 800 GeV) and NA50 [3] (E p = 450 GcV). As 
already seen (e.g. in Fig. 7), the measured pp de- 
pendence of the nuclear modification effects is very 
similar for HERA-B and E866. In Fig. 11 the E866 
and HERA-B measurements are seen to be compati- 
ble within statistical and systematic uncertainties in 
the overlap region. The NA50 results are based on 
lower energy collisions and are systematically below 
both HERA-B and E866. At lower values of x F , the 
HERA-B o.(xf) measurement indicates a reversal of 
the suppression trend seen at high xp: the strong 
suppression established by previous measurements at 
high xp turns into a slight tendency towards enhance- 
ment in the negative xp region. 

The dependence of J/ip production in hadron- 
nuclcus interactions on xp has been modeled by 
Vogt [4, 5]. Nuclear effects caused by final-state ab- 
sorption, interactions with co-movers, shadowing of 
parton distributions, energy loss and intrinsic charm 
quark components are described separately and in- 
tegrated into the model. It is further assumed that 
the cc pair is subject to more severe energy losses if 
produced in a color octet state. Four curves from 
this model which differ in their descriptions of nu- 
clear Parton Density Functions (nPDF) and energy 
loss are shown in Fig. 1 1 . All calculations shown here 
were done for the center-of-mass energy of HERA-B 
at Vs = 41.6 GcV. The nPDF distributions of Eskola, 



a 1.1 



t HERA-B 2008 e + \i 
□ E866 2000 




X F 



Figure 11: Measurements of a as a function of xf by HERA-B (filled triangles, plotted with total and 
statistical uncertainties), E866 (y/s = 38.8 GeV) [2] (empty squares) and NA50 (y/s = 29.0 GeV) [3] (empty 
triangles). The curves were calculated by Vogt [4, 5] based on three different nuclear parton distribution 
functions: (EPS [37], EKS [33, 39] and HKN [38]) and two models of initial state energy-loss: GM [7] and 
BH [8]. For all approaches, energy loss, intrinsic charm and shadowing are taken into account. 



Kolhincn and Salgado (EKS) [33] describe the scale 
dependence of the ratios of nPDFs of a proton inside a 
nucleus to those of a free proton within the framework 
of lowest order leading-twist DGLAP evolution [34] 
by evolving the initial PDF from the CTEQ4L [35] 
and leading order GRV [36] parameterizations. An 
improved leading-order DGLAP analysis of nPDFs 
including next to leading order calculations has been 
published recently by Eskola, Paukkunen and Salgado 
(EPS) [37]. In another approach by Hirai, Kumano 
and Nagai (HKN) [38] , nuclear structure function ra- 
tios F2/F2 and Drell-Yan cross section ratios are 
analyzed to obtain nPDFs. The HKN analysis shows 
weak anti-shadowing at negative xf- 

Initial state energy loss as described by Gavin and 
Milana (GM) [7] and modified by Brodsky and Hoyer 
(BH) [8] is based on a multiple scattering approach 
that essentially depletes the projectile parton momen- 
tum fraction as the parton moves through the nucleus. 
Both quarks and gluons can scatter elastically and 
therefore lose energy prior to the hard process result- 
ing in an effective reduction of J/ip production for 



x F > 0. 

The measurement of HERA-B shows that a in- 
creases with decreasing xp and suggests enhanced 
J ftp production for x F < -0.1. The HERA-B data 
favors the nPDFs of EPS and HKN over EKS. The BH 
description of energy loss is clearly ruled out. None 
of the variants of the Vogt model give a satisfactory 
description of both HERA-B and E866 data. For ex- 
ample, while the HKN curve is compatible with most 
of the HERA-B data points at negative xf 1 it lies 
significantly above the E866 points and furthermore 
fails to adequately describe RHIC data [37]. 

Another theoretical model by Borcskov and 
Kaidalov [6] , formulated in the framework of reggeon 
phenomenology, predicts an anti-screening effect in 
the region of negative xp- An important ingredient 
of their model is the assumption that a colorless state 
containing c and c quarks which has some probability 
of projecting into a charmonium state is produced and 
propagates through the nucleus. The colorless state 
is of large size, possibly consisting of DD or D*D* 



Vt (GeV/c) 


a 
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a a 
0.0 


0.0 


0.906 ± 0.006 ± 0.021 
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1 099 ± 026 ± 028 


/I R 
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0.4 


1 139 ± 056 ± 038 


Xp 




a 


min 


max 




-0.34 


-0.26 


1.036 ± 0.034 ± 0.042 


-0.26 


-0.22 


1.012 ± 0.023 ± 0.030 


-0.22 


-0.18 


1.031 ± 0.014 ± 0.023 


-0.18 


-0.14 


1 A1C_I_AA1A_I_A A9n 

1.U10 It U.U1U It u.uzu 


-0.14 


-0.10 


0.994 ± 0.008 ± 0.018 


-0.10 


-0.06 


0.978 ± 0.006 ± 0.018 


-0.06 


-0.02 


0.967 ± 0.005 ±0.017 


-0.02 


0.02 


0.967 ± 0.006 ± 0.018 


0.02 


0.06 


0.962 ± 0.008 ± 0.021 


0.06 


0.14 


0.947 ± 0.015 ± 0.028 



Table 6: The parameter a as a function of pr and xf- 
Statistical and systematic uncertainties are indicated 
separately. They were extracted from the measure- 
ment using Eq.s 13 and 17. 



mesons, and therefore has a large interaction cross 
section. As it propagates through the nucleus it inter- 
acts and loses energy. Ultimately, the observed J/ip 
mesons are projected out of the energy-depleted col- 
orless state. The measurements of HERA-B are qual- 
itatively compatible with the calculations described 
in BK [6]. 



6 Conclusions 

HERA-B has performed the first determination of the 
nuclear dependence of J/ip production kinematics at 
negative xf in proton- nucleus collisions. The ana- 
lyzed data samples were obtained in collisions of pro- 
tons from the 920 GeV HERA-proton beam with car- 
bon, tungsten and titanium targets. The J/ip mesons 
are observed in both dimuon and dielectron decay 
channels. The comparison of results from the two 
channels affords some additional control over system- 
atic uncertainties arising from triggering and selection 
procedures. 

The measurement covers the kinematic range 
-0.34 < x F < 0.14 and p T < 5.4GeV/c. The mea- 
sured dN/dpr distribution is seen to become broader 



with increasing atomic mass number as has already 
been observed by experiments at lower center-of-mass 
energies [1, 3]. The data indicates that the dN/dxF 
distribution also tends to become broader and that 
its center moves towards negative xf values with in- 
creasing A. 

The dependences of the nuclear suppression pa- 
rameter, a, on pt and xf are also presented. The 
a parameter is seen to increase with increasing pt 
in agreement with data from E866 [2]. In the xf 
region of overlap of the two experiments, the two a 
measurements are mutually consistent. As if de- 
creases, a increases and becomes greater than 1 be- 
low xf ~ —0.15, although the data remains compat- 
ible with a value of 1 to within 2 a. Thus instead 
of the strong suppression observed at high positive 
xf, HERA-B measured no suppression or a possi- 
ble enhancement of J/ip production at negative xp- 
Hard-scattering based models [4, 5] have difficulty 
simultaneously accommodating the HERA-B, E866 
and RHIC measurements, while the reggeon-inspircd 
model of Boreskov and Kaidalov [6] is in qualitative 
agreement with the data. 
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